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Abstract 

We use the Renormahzation Group method to study the Bose-Einstein condensation of the 
interacting dilute magnons which appears in three dimensional spin systems in magnetic field. 
The obtained temperature dependence of the critical field Hc{T) — Hc{0) ~ is different from 
the recent self-consistent Hartree-Fock- Popov calculations (c ond-matt 0405422|) in which a T'^/'^ 
dependence was reported . The origin of this difference is discussed in the framework of quantum 
critical phenomena. 
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I. INTRODUCTION 



Recently, the quantum spin system in magnetic field became of great interest because 
the experimental data on the ladder systems showed the possibility of a quantum phase 
transition (QPT) driven by the magnetic fieldQ]. The QPT in the spin systems have been 
treated using the quantum rotor modelP]. However, the transition driven by a magnetic 
field is difficult to be described in such a formalism. 

The occurrence of the magnetic order in the spin-gap magnetic compounds has been 
interpreted as a Bose-Einstein condensation (BEC) of magnons. At the present time it 
is known that the BEC has been discovered in ultra-cooled dilute atomic gases, but these 
experiments present various limitations. The analogy between a quantum spin system which 
presents long-range order and an interacting Bose gas which presents BEC is well known 
for a long timej^. In the case of quantum spin system is possible to tune the density of 



magnons by magnetic field 
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or pressure to observe BEC in these systems. 



In this paper we formulate the theory of magnon condensation using the Renormalization 
Group method (RNG) applied for bosonic systems J, 0, ^ Q, 0] and recently reconsidered 
for the spin systems|lO|]. 

The outline of this paper is as follows. In Sect. II we present the basic experimental 
evidence obtained on the material class XCuCl^ {X = Tl, K, NH4) and the main results 
obtained with the mean field theory. The model based on these data will be presented in 
Sect. III. The Sect. IV is devoted to the RNG formulation using results from Ref. ^ 
in order to calculate the low temperature quantum critical properties of the system. Sec. 
V will present the main results of the paper, the relevant thermodynamic quantities next 
to the quantum critical point. Finally, in the last section we will discuss the results and 
compare with other theoretical approaches and experimental evidence. 



II. EXPERIMENTAL EVIDENCES 



In this section we present the basic experimental data obtained on the XCuCl^ com- 
pounds which gave us the possibility to elaborate a model for the critical behavior of this 
QPT clS cl BEC driven by the magnetic field. 

The compound which offered the possibility of performing many measurements is 
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TlCuCl-i which is composed of a chemical double chain Cu2CIq. The magnetic suscep- 
tibility, measured as function of temperature T for three different directions of the magnetic 
field, exhibit a broad maxima at T = 38/^, decreasing to zero with the temperature de- 
creasing 11]. This result indicates that the compound has an excitation gap ^ = l.hK 
above the singlet ground state. This gap may by attributed to the antiferromagnetic dimer 



coupling in the double chain [1^ . The magnon dispersion has been investigated theoretically 
jlS^ and experimentally by inelastic neutron scattering Q|. The experiments showed that 
the magnon modes are split into three by magnetic field with the splittings proportional to 
the field , and the lower modes becomes soft at the critical field He. In a magnetic field H 
the chemical potential is A — gfisH (where g = 2) and He = 

An important result has been obtained by Sherman et al. by sound attenuation in 
TlCuCl^ in magnetic field at low temperatures. The occurence of the sharp peak in the 
sound attenuation near the BEC critical temperature and the Drude form of the sound 
dumping suggested a constant weak magnon- magnon coupling. 

The sound attenuation suggested a magnons dispersion like tu^ = -|- ^^j^ (here J is 
the exchange interaction) for small wave vector " /c" , approximation valid if Tc and A are of 
the order of fc^Tc = 7K and fc^A = 7.5 as for the case of TlCuCl^. 

The existence of a induced magnetic field-spin transition has been proposed by Giamarchi 

n 

and Tsvelik 16] considering as a possible mechanism the BEC of the soft mode. The Hartree- 
Fock approximation using the Hamiltonian from the theory of BEC in dilute bosonic gases 
with effective chemical potential = iXBg{H — H^) has been applied |l7| to calculate the 
temperature dependence of magnetization. Near the critical temperature Tc this approxi- 
mation breaks down and the magnetization is expected to behave like M ~ (T — T^)^ where 
/9 = |. The temperature dependence of the critical field has the form [Hc{T) — Hc{0)] ~ T* 
where $ = | but the best fit is obtained for $ = 2.2. 

The system has some important characteristics. First, we menton the small three- 
dimensional (3D) character of the system which leads to a small gap. Another important 
feature is the decreasing the number of thermal excitations with the increasing of the mag- 
netic field, and as a consequence the driving the condensate to T = behavior. 

The form of the magnon dispersion (called "relativistic") changes the dynamics of the 
system and we expect $ = 2. This conjecture is on of the most important result of Ref. Q| 
and we will use it for the model of the magnon condensation. 
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III. MICROSCOPIC MODEL 



We introduce the description of the magnon condensation in magnetic field using the 
action: 



where: 



^eff - Siff + Siff (1) 



Si% = ^ E X-\k)\m\' (2) 



Si% = ^ E • • • E 'Pik,) . . . <Pih)5ih + + (3) 



-,(4) _ Uo 

ki kn 

Here we introduced the notations k = (k, u;„), Un being the bosonic Matsubara frequency 
and 

In Eq.(I2I) x{k) is the magnon propagator and uq the bare couphng constant. Using [1^ 
we write the magnon propagator as 

+ + ro 

where tq = A — ^bqH and H is the external magnetic field. 



IV. RENORMALIZATION GROUP EQUATIONS 

The model has the dynamical critical exponent z = 1 and for this model d = 3. The 
Renormalization group equations in one- loop approximation at finite temperature T for 
d = 3 have the form 



- m (5) 



du{l) _ {n + 8)K3 2 
dl ~ 8 



dl 

'''^\\i)FMi),m] (6) 



^-M/)+(!i±|l^.(OW),T(0] (7) 

where n is the number of components of the fiuctuation field 7^^3=2^. We denote 

by u{l = 0) = tio,and r{l = 0) = tq and the functions Fi 2 are characteristic functions for 
the model We approximate F2[r{l),T{l)] ~ ^ in the low temperatures limit and the 
solution of Eq. has the form: 



where Cq = ^"''^^^^^ and K^, = The physics near the QCP is described by the scahng 
field tr{l) defined as: 

T7 -I- 2 

^.(0 = ril) + -j^KMl) (9) 
The expression of tr{l) has been calculated in Ref. [3| for d=3 as: 

4 "'^i '^■^ eVn-)-l 



tr{l) = e'^^W k{0) + /' dx I (10) 



where Ar{l) has the expression: 

T7 -I- 2 / 

A.(/) = 2/-^ln(f + 1) (11) 

Using these results we will calculate the thermodynamic quantities in near the critical point. 
The basic idea of the method is to stop the renormalization procedure close enough to this 
point that the system can see the infiuence of the quantum effects. This matching condition 
is in fact equivalent to the stoping of the renormalization at the scale I = Z*, where /* ^ 1. 
This value is obtained from the condition [l^: 

uir) = 1 (12) 

In the approximation I* > 1 and T(/*) > 1 we obtain from Eq.()ll|) the relation: 

exp(r) = [t,„(T)]ir (13) 
where is the eigenvalue of the relevant parameter r and is given by the expression: 

A^ = 2,e = 0(rf = 3) (14) 
and troiT) is given by the relation jl^ : 

tro{T) = ro - roc + -^^'^o'^^ (15) 



Following the method from [3, Q we calculate /* as: 



/* = - In (16) 
2 T ^ ^ 
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V. THERMODYNAMIC QUANTITIES 



Using these results we will calculate the relevant thermodynamic quantities near the 
quantum critical point, but in the disordered state. First, we define the critical line by 
t{T) = and we get: 

H,{T) - H,{0) = CoT^ (17) 



where Cq oc uq and He = The temperature dependence of the number of magnons is 

n{T) = exp{-3n J -^M^in] (18) 



defined as: 

(2vr) 

This equation gives for n{T) a temperature dependence of the form: 



n(T) oc (19) 

This result is in agreement with the behavior of magnetization at very low temperatures. 

From Eq. (jl5j) we calculate the critical line in the (ro,T) plane using the condition 
troiT) = 0. This gives: 

roc{T) = roc - ^^^^«oT' (20) 



and for tq < tqc we get the general equation 



Te(ro) 



128 



[n + 2)Uo 

Using now the definition (see Ref. Q|) of roc 



1/2 



{roc -To)'/' (21) 



in + 2) 

we obtain 

Tc{H) \H - Hc\'/" (23) 



with a = 2 and 

(n + 2) 
327r 



The specific heat can be also calculated for the important case vq ^ tqc but T T+(ro) 
using the singular part of the free energy reported in [lO] as: 

Fs{T) ~ [T-Tc{H)]^ \ \n{T-n{H))\[''^] (24) 
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and it give us a logarithmic behavior: 

^^|ln(T-T.(/7))|HSi] (25) 

We will discuss these results in connection with the existent theoretical approaches and the 
experimental data for TlCuCl^. 

VI. DISCUSSIONS 

The occurence of BEC induced by the magnetic field in a spin-gap system has been 
predicted by Giamarki and Tsvelik (1(tI | in order to explain the three-dimensional ordering in 
coupled ladders. Nikumi et al. [17] applied Popov theory (see for example Q]) for the BEC 
obtaining a temperature dependence for magnetization in agreement with the experimental 
data, of the form T^/^. 

Howeve, a. i. wa. mentioned i„ Ref. Q .he cn.ical expose. * of the cnt.cal field 
[Hc{T) — -f^c(O)] ~ T* was obtained as $ = 3/2, but the experimental data shows $ = 2.2. 

Recently Misguich and Oshikawa reconsidered the self-consistent Hartree-Fock-Popov 
(HPF) method using a realistic dispersion for magnons jisl and they reobtained $ = 3/2. 
They also calculated the specific heat and obtained a A-shape. 

We would like to mention several basic points related to this problem: 

• in the HPF the effect of fluctuations have been neglected and it appears normal that 
thermodynamic quantities on the critical line have wrong critical exponents. 

• the changing of the dispersion law for the magnon can improve the HEP results, but 
does not affect the basic approach. 

• both papers and j3l considered the QPT as driven by the temperature and in 
fact it is induced by the magnetic field. 

These observations justify the application of the RNG method in the version proposed by 

n 

Caramico et al. to study the magnon condensation. Our main result can be summarized 
as: 

• we obtained the critical exponent $ = 2, which is in a better agreement with experi- 
mental data ($ = 2.2). 



• the magnetization calculated in our approach is also T^/^ dependent. 

• the critical temperature Tc{H) has a dependence given by Eq. (j^BJ with a = 2 in 
agreement with The numerical results have been analyzed in j2dl |. 

• the specific heat calculated by RNG method has a logarithmic correction as was ex- 
pected. The T dependence of is not relevant if we take the approach of a magnetic 
field driven OPT, but we have it to the dependence of H(T). Such an experiment has 
been done and we can notice a small anomaly in this dependence. But a final 
decision about the H{T) dependence of our Eq.(25) cannot be taken. 

The agreement between the experimental data and the theory presented in appears in 
our opinion due to the fact that the authors considered, on a very narrow interval a linear 
dependence between the critical field Hc{T) and the critical densitync. This approximation 
can be valid for ric < 0.002 which describes a part of the real temperature dependence of 

VII. ACKNOWLEDGMENTS 

One of the authors (M. C.) thanks Luigi De Cesare for numerous correspondences helping 
him to understand the RG methods developed by the Salerno group(A. Caramico D'Auria, 
L. De Cesare and I. Rabuffo) for this problem. 



[1] S. Sachdev, Science 288, 475 (2000). 

[2] S. Sachdev, Quantum Phase transitions, (Cambridge Univ. Press, 1999). 

[3] T. Matsubara and H. Matsuda Prog. Theor. Phys. 16 569(1956). 

[4] D. S. Fisher and P. C. Hohenberg Phy. Rev. B 37 4936 (1988). 

[5] A. I. Uzumov, Introduction to the Theory of Critical Phenomena (World Scien- 
tific, Singapore, 1993). 

[6] S. Sachdev, T. Sentil,and R. Sankar Phys. Rev. B 50, 258 (1994). 

[7] A. Caramico DAuria, L. De Cesare, and I. Rabufo Physica A152 (1997). 

[8] M. Crisan, D. Bodea, I. Grosu and I. Tifrea J. Phys. A; Math. Gen.35, 239 (2002). 



8 



A. A. J. Schakel, Icond-mat ()3()1()5()1 

A. Caramico D'Auria, L. De Cesare, and I. Rabufo Phys. Lett. A308, 208(2003); Physica 
A327,442(2003). 

K. Takasu, W. Shimramura and H. Tanaka, J. Phys. Soc. Japan. 66, 1611 (1997). 

A. Oosawa, M. Ishi and H. Tanaka J. Phys; Condens. Matter 11 265, (1999). 

M. Matsumoto, B. Normand, T. M. Rice and M. Sigrist Phys. Rev. Lett. 89, 077203 (2003). 

N. Cavadini, G. Heigold, W. Henggler, A. Furrer, H-U Gudel, K. Kramer and H. Mutka Phys. 

Rev. B63, 172414 (2001) 

E. Ya. Sherman, P. Lemmens, B. Busse, A. Oosawa, and H. Tanaka Phys. Rev. Lett. 91, 
057201 (2003). 

T. Giamarchi and A. M. Tsvehk Phys. Rev. B59, 11398 (1999). 

T. Nikuni, M. Oshikawa, A. Oosawa and H. Tanaka Phys. Rev. Lett.84, 5868 (2000). 

H. Si and A. Griffin Phys. Rep. 304, 1 (1978) 

G. Misguich and M. Oshikawa cond-mat /0040542p 



O. Nohamadani S. Wessel B. Normand and S. Haas , cond-mat / 030712'6| ; B. Normand, M. 



Matsumoto, O. Nohadani, S. Weesel, S. Haas, T. M. Rice, and M. Sigrist J. Phy: Condens 
Matter, 16, S867 (2004) 



9 



